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ABSTRACT 


In recent decades following a lot of droughts, many changes have been made in the quantity and quality of 
the country’s water resources. This factor has caused many uncertainties in the management of the country’s 
water resources. The purpose of this study is to improve the understanding of drought effects on the quantity 
and quality of water resources in Lorestan province from the years 2008 to 2018 by coherence and cross 
wavelet method. To achieve this goal, first, drought assessment according to precipitation data has been 
examined using (SPI) index, and then the effect of drought on Khorram river runoff is analyzed. Then, the 
global index of water and the drought impact on this index in the Khorram River were evaluated. Also, the 
global index of water quality (WQI) and the impact of drought on this index in the Khorram River are evaluated. 
The results of coherence and cross wavelet indicated the relative effect of precipitation with a wavelet 
coherence coefficient of 0.6 on changes in water runoff in the Khorram River is of degree first importance. 
Also, the relative impact of drought with a wavelet coherence coefficient of 0.4 changes in water quality of 
Khorram River has been more than other factors. Therefore, climatic factors in reducing the water runoff of 
the Khorram River from factors other are more important. Also, the results showed human factors as the most 
important ones in water quality changes in the Khorramabad River. 
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1. INTRODUCTION 


ivers are a small part of the world’s running water 
R« carry about 32 to 37 cubic kilometers of water 

to the ocean each year. Rivers are one of the main 
sources of water supply for various purposes such as 
agriculture, drinking, and industry. Iran is located on the 
world’s drought belt by being in the position of 25 to 40- 
degree north latitude and for a long time and now is facing 


the problem of drought. This subject, besides the 


occurrence of severe and prolonged droughts, created 
conditions for the country, which has faced serious 
challenges. Numerous studies have been conducted to 
identify the effects of drought on the quantity and quality 
of water resources. For example, Researchers studied the 
statistical analysis relating variations in groundwater 
levels to droughts on Jeju Island, Korea. This study 
indicated that PCA could be a powerful tool for 
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summarizing large datasets to select index wells that are 
vulnerable to droughts, which would significantly reduce 
the expense of monitoring programs (Jung et al. 2021) [1]. 
In this regard, Komasi et al. (2016) examined the routing 
and classification of factors affecting the reduction of 
groundwater levels by using coherence and cross wavelets 
in the Silakhor plain. The result of this study showed that 
the runoff time series, which represents human effects, has 
a greater effect on reducing the aquifer level of Silakhor 
plain than the precipitation and temperature time series, 
which represents climate change [2]. Also, Soleimani 
Motlagh et al. (2015), in a study on the aquifer resources 
of Lorestan province, showed that drought has a 
significant effect on water quality variables such as Ci, Na, 
TDS, and EC. One of the important results of this research 
is that the water quality based on Schuler and Wilcox 
diagrams has changed due to drought, and the quality of 
drinking water has decreased [3]. Wolff et al. (2021) 
studied the impact of the 2018 drought on pharmaceutical 
concentrations and the general water quality of the Rhine 
and Meuse rivers. The results showed that low flow 
combined with high temperatures resulted in a general 
deterioration of surface water quality of both the Meuse 
and Rhine rivers during the 2018 drought. Also, the trend 
of declining groundwater levels continues to decline, and 
the continuation of this trend will further reduce the runoff 
of wells, canals, and springs [4]. Han et al. (2020) studied 
the effects of vegetation restoration on groundwater 
drought in the Loess Plateau, China. Results indicated that 
the normalized difference Vegetation Index (NDVI) 
during 2003 — 2015 in the LP was growing rapidly; 
meanwhile, groundwater storage significantly decreased 
(p < 0.01), and groundwater drought intensified in terms 
of its area and intensity [5]. Piravi et al. (2015) a study 
evaluated and analyzed the effects of drought on the water 
quality of rural wells in the Kashfar area of Mashhad plain. 
The results showed that the groundwater level decreased 
quantitatively and qualitatively from 2007 to 2011, and a 
significant relationship exists between changes in water 
level and water hardness [6]. The results of a study 
examining the groundwater drought index in the United 
States have shown that the SPI drought index with a delay 
of 12 and 24 months has the highest correlation with the 
SWI index (Tatiana 2019) [7]. Zhao et al. (2021) 
investigated quantitative analysis of nonlinear climate 
change impact on drought based on the standardized 
precipitation and evapotranspiration index. The results 
indicated that drought on the SNP has been mitigated 
marginally during 1961-2016, mainly during spring, 
winter, the growing season, and on the annual timescale 
[8]. In another study, neural network and regression 
models were used to predict drought in the Polish- 
Ethiopian river basin. The results showed that the linear 
regression model is better than the artificial neural network 
(Mishra et al. 2012) [9]. Kim et al. (2019) investigated the 
quantitative vulnerability assessment of water quality to 


extreme drought in a changing climate. The results of this 
study showed that due to severe droughts, the water quality 
of the Nakdong River has decreased [10]. Ghaffari et al. 
(2021) studied the Spatial and temporal variation of 
groundwater quality around a volcanic mountain in 
northwest Iran. The results indicated that groundwater 
quality in the studied area is deteriorating [11]. On the 
other hand, Nourani et al. (2015) investigated the changes 
in hydrological processes using the entropy wavelet 
criterion in Lake Urmia. Using this criterion, they 
introduced the reduction of water flow time series 
fluctuations in the Urmia Lake basin as the most important 
factor in reducing the water level of this lake [12]. To 
investigate the complexity of the entropy method, he used 
several scales for rainfall and runoff time series and 
observed that the results obtained at higher time scales 
were different from the results obtained at lower time 
scales (Chou 2014) [13]. Feng et al. (2022) investigated 
climate change impacts on concurrences of hydrological 
droughts and high-temperature extremes in a semi-arid 
river basin of China. Results showed that the frequency of 
CHDHEs increased by 160% from 1961 to 1988 to the 
recent period 1989-2016 [14]. Evan & Timothy. (2021) 
studied the effects of climate and land use changes on the 
water quantity and quality of coastal watersheds of 
Narragansett Bay. Results showed significant effects of 
climate change and land-use change on the watershed, 
with demonstrated impacts on sediment loading, organic 
N, organic P, and nitrates. Climate impacts were much 
more significant than land-use effects, but land-use 
impacts displayed greater regional variation [15]. Mohd & 
Abhishek. (2020), the investigated groundwater quality 
assessment in the Lower Ganga Basin using entropy 
information theory and GIS. The results obtained in this 
paper are useful for identifying suitable sites for tube wells 
and bore-wells for drinking and agricultural purposes and 
also for developing an effective strategy to avoid further 
contamination of groundwater aquifers in the region [16]. 
Kubicz et al. (2021) studied the effects of drought on 
health risks 
contamination. This study found that in some cases, the 


environmental posed by groundwater 
occurrence of drought did not cause an increase in the non- 
cancerogenic threat expressed by the hazard index [17]. 
Hassanzadeh et al. (2021) investigated the effect of climate 
change on the decline of Khorramabad groundwater level 
using an entropy wavelet. The results showed that climatic 
factors (precipitation) had a greater effect on reducing the 
[18]. 
background, the analysis of data and time series by wavelet 


groundwater level According to the research 
transform with high speed and accuracy was performed, 
and unlike other methods of signal analysis, the trend of 
sudden changes, breakpoints, and discontinuities can be 
identified. Due to the comprehensiveness and suitability of 
the SPI index for drought monitoring, therefore, from SPI 
standard precipitation index to drought monitoring and 
also, in order to analyze and evaluate the quantity and 


88 


J. Civil Eng. Mater.App. 2022 (June); 6(2): 87-98 


quality of Khorram river water with different factors have 
been used, cross and coherence wavelet and the global 


2. MATERIALS AND METHODS 
2.1. CASE STUDY 


Khorram River is in Khorramabad city in Iran, which is 
located between the lengths of 47 °55’ to 48 °50’ east and 
latitudes 40 °32’ to 34 °20’ north and at an altitude of 1147 
meters above sea level. Also, the length of this river is 
about 110km, of which about 5km passes through the city 
of Khorramabad. This river joins the Kashkan River in the 
place of two waters and after a long route and joining rivers 
such as Cholhol and Madian River in the place of 
Gaumishan Bridge; it is connected to Seymareh River. The 


quality index WQI. 


city has a temperate and semi-moderate Mediterranean 
climate, with favorable precipitation, especially in spring. 
The water of the Khorram River enters through four 
stations Bahramjo, Chamanjir, Doabvisian, and Karganeh. 
The highest and lowest runoff of the Khorramabad River 
in the years 2008 to 2019 was in 2018 and 2012, 
respectively. Table 1 shows the statistical characteristics 
of the study area. 


Table 1. Statistical information of the study area from 2008 to 2018 


Parameter Minimum 
Precipitation (mm) 0 
Temperature (degree Celsius) -l 
Evaporation (mm in month) 0 
Runoff (Cubic meters per second) 146.5 


2.2. STANDARD PRECIPITATION INDEX (SPI) 
The SPI index is a powerful tool to process the 
precipitation data, and its purpose is to assign a numerical 
value to precipitation through which areas with different 
climates can be compared together. One of the advantages 
of SPI is that the calculation of SPI is based on rainfall data 
and does not depend on soil moisture conditions; another 
advantage is that this index is not affected by topography 


Maximum Average Standard Skewness 
deviation 
1255 38.1 46.5 IES) 
40.3 19.2 8.5, 0.8 
360.5 i) 140.3 0.2 
452.1 238.2 144.7 2.2 


(Komasi et al. 2013) [19]. In general, this index is defined 
to express drought as follows: when the SPI is permanently 
negative and reaches a value -1 or less and when the value 
is positive. Therefore, positive values indicate above- 
average rainfall, and negative values indicate less rainfall 
than average rainfall, which is shown in Table (2) (Komasi 
et al. 2013) [19]. 


Table 2. Drought classification based on SPI index (Hassanzadeh et al. 2012) [20 


Drought intensity 


extreme wet 


very wet 


Moderate wet 


Near to normal 


Moderate drought 


very drought 


extreme drought 


The main data of the SPI index are rainfall data of rain 
gauge stations. After ensuring the homogeneity and 
randomness of the monthly data, the time series is formed 
in the intervals of 6 and 12 months, and its time series is 


Index values 


2 and more 
1.5 to 1.99 
1 to 1.49 
-0.99 to 0.99 
-1.49 to -1 
-1.99 to -1.5 
-2 and lees 
fitted with the gamma distribution, the probability or 


plenty function of which is given in Equation | (Mckee et 
al. 1993) [21]. 
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Xo- 1 e*/B (1) 


Bree a) 


In this relation x=0 is the amount of precipitation, o> 0 is (a) is the gamma function represented by Equation 2. 
the shape parameter, B> 0 is the scale parameter and T 


(a) = fp yo te dy (2) 
In equation (2), the parameters a and f related to the gamma estimated a and B coefficients using the optimal maximum 
density function are estimated for each station and for each proofing base on equations (3) to (5). (Mckee et al. 1993) 
time scale, and for each month of the year. Mckee et al. Aue 

1 4A 
a= 7A (1 + f1+ “) (3) 
= In(X 

a =ingy — 22900 a 

B=- (5) 
In the above equation, n is the number of precipitation station; this probability can be converted to an incomplete 
observations, and x is the average precipitation for several gamma function according to equation 6 with the 
months. The calculated parameters to find the cumulative assumption t=x/B. 


probability of precipitation used for a specific time in each 


G(X) J g(X)dX = a fo te-Y eat (6) 


Since the gamma function is not defined for x=0, the probability is calculated as equation 7. 
rainfall distribution has a value of zero the cumulative 


H(X) = q+ (1-q)G(X) (7) 
Q is the probability of rainfall zero, and m is the number the product of m divided by the total number of data (n), 
of zeros in the rainfall time series, which estimates (q) as then by having H (X) and relations 8 to 11, SPI is obtained. 
Coty t+Cot 
onl = -|t- al 0<H(X) < 0.5 (8) 
_ CotCyt+Cot 
SPI = +|t— el 0.5 < H(X) <1 (9) 
In (=) 0<H(X) <05 (10) 
= + /In — am) 0.5 < H(X) <1 (11) 
In these relation, the coefficients are constant in relationships 8 to 11 (Mckee et al. 1993) [21]. 


Co,C1,C2,d1,d2,d3_ which should be placed from table 3 


Table 3. Coefficients values of SPI calculation formulas (Mckee et al. 1993) [21 
Coefficient dy dy ds Co Cc, C, 
Value 1.43 0.18 0.01 2.51 0.80 0.10 


2.3. CROSS WAVELET TRANSFORMATION 


The main purpose of a cross wavelet is to obtain a diagrams identify periods that provide areas with high 
complete time-frequency representation of a local and examined wavelet spectra examined. According to 
temporary event that varies on time scales. Cross wavelet equation 12 with any desired mother wavelet, like Morlet’s 
transformation diagrams are examined to identify periods mother wavelet, equation 13 can estimate the wavelet 


that offer regions with high wavelet spectra. Cross wavelet 
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transform for the time series of each hydrological data x(t) 
(Labat. 2010) [22]. 


C¥ (ab) = fxOQZyE) 


Wo(n) = M~*/4e!@one-n’/2 


Yo is a function of the mother wavelet, e exponential 
function and w frequency without dimension n time are 
dimensionless, and also mark "=" refers to a mixed mating 


of the mother wavelet. The parameter "a" is expressed as 
a scale factor if it is a >1, the time series expands along the 
time axis, and if a<1 the time series contracts along the 
time axis. Parameter "b" is used as a position factor and 


W;,(a,b)= C ‘ (a, b)C ‘i (a, b) =|cx(a,b)|? 


This wavelet spectrum can be averaged in time, which is 
generally defined as the average wavelet power spectrum 
and allows to specify the scale specification (Torrence & 
Compo 1998) [23]. The vacillation alternation period 
specification is determined using the overall wavelet 


Wxy(a,b)= Cy (a,b) CY (a,b) 


Which Cy (a,b) and C*y (a,b) are the continuous-time 


series wavelet coefficients x(t) and mixed, the wavelet 


2.4. COHERENCE WAVELET TRANSFORMATION 


The coherence wavelet spectrum was inappropriate to 
express the interrelationship between the two processes, 
and hence the use of coherence wavelet transforms in time 
series processing which is more appropriate to find a 
criterion of correlation between two-time series at 
different frequencies and time (Nourani et al. 2015) [24]. 


(12) 


(13) 


allows you to study the time series x(t) around time b. the 
concept of wavelet transform can be used to investigate the 
relationship between two different time series related to 
two separate hydrological processes. For this purpose, the 
wavelet spectrum W,(a,b) of the time series x(t) is similar 
to Fourier analysis and is defined by the absolute value of 
the wavelet coefficient. 


(14) 


spectrum. Similar to the Fourier coherence spectrum and 
wavelet coherence spectrum W,,(a,b) is defined between 
two different hydrological time series x(t) and y(t), as 
follows. 


(5) 


coefficient is y(t). The wavelet spectrum averaging 
technique is used to express the mutual covariance of time 
series x(t) and y(t) and their distribution at different scales. 


In this regard (Torrence & Webster 1999) [25] suggested 
that the wavelet communication be determined using 
smooth wavelet spectrum estimation. Smooth wavelet 
spectrum SWxx(a,b) and wavelet coherence spectrum 
SWxy(a,b) are defined as below: 


SWxx(a, b) = (aes W*x(a,b) Wx(a,b)dadb (16) 
t+6/2 + 
SW,y(a, b) = (hare W'x(a,b) Wyy(a,b) (17) 


In this regard, 5 represent the size of the two-dimensional 
filter (Luterbacher et al. 2002) [26]. Finally, the criterion 


WC (a,b) _ |SWxy(a,b)| 


2.5. WATER QUALITY INDEX (WOQI) 

The water quality index WQIT is one of the techniques used 
to evaluate water quality. This method was first proposed 
by Horton in the year 1970. This factor is usually obtained 
with the number of general parameters of water, including 
dissolved oxygen, acidity, hardness, water-soluble solids, 
temperature, turbidity, and the electrical conductivity of 


VlSwxx(a,b)|.|SWyy(a,b)| 


of coherence wavelet transform relation can also be 
defined similarly to Fourier coherence as follows: 


(18) 


water, nitrite, nitrate, chlorine, and some of the main ions. 
Some studies use this statistical technique to evaluate the 
water quality index by using the weight score of each 
parameter analyzed (Toledo et al. 2002) (table 4) [27]. In 
this classification, waters with WQI less than 50 are in the 
very good category, 50 to 100 are in a good category, 100 
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to 200 are in the poor category, 200 to 300 are in the very 
poor category, and more than 300, are in the class of 
unsuitable drinking water. In this study, the parameters 
PH, TDS, Ca, Mg, Na, K, HCO3, Cl, SO4, TH, and EC 
were used for the WQI index. In WQI index calculation, 


Where, w; weight of each parameter n number of parameters. 


the first step is weighting for each parameter. According to 
its relative importance, specific weight is assigned (table 
4). The second step is the calculation of relative weight 
based on equation (19): 


(19) 


Table 4. Weight ratio of chemical parameters (WHO Standard) 


Chemical parameters Unit 
PH - 
TH Mg/l 
EC p. mhos/cem 
Ca Mg/l 
Mg Mg/l 

K Mg/l 
Na Mg/l 
HCO3 Mg/l 
Cl Mg/l 
S04 Mg/l 
TDS Mg/l 
Total - 


The third step is to calculate the quality rate scale. This 
scale (qi) is calculated by dividing the concentration of 


a= © 100 


Which in Ci concentration and Si chemical standard each 
parameter in the water sample is in Mg/I. then to estimate 


SHi=W gi 


WQI= 2, Sh, 


WHO Standard Weight Weight ratio 

6.5-8.5 4 0.129 
200 2 0.064 
250 3 0.096 
75 2 0.064 
50 1 0.032 
12 2 0.064 
200 2 0.064 
120 3 0.096 
250 3 0.096 
250 4 0.129 
600 5 0.161 


31 


each parameter in each water sample by the standard value 
of that parameter according to equation (20). 


(20) 


WQI, Si for each chemical parameter is determined by 
equations (21) and (22): 


(21) 


(22) 


Table 5. Water quality classification based on WQI index (Ramakrishanaiah et al. 2009) [28] 


Quality class 
Unsuitable 
Very weak 

Weak 


good 


Excellent 
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Value of the index obtained 


300 


200-300 


100-200 


50-100 
<50 
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3. RESULT AND DISCUSSION 
3.1. DROUGHT INDEX RESULT 


To calculate the drought index, from precipitation data 
obtained, the meteorological organization of Lorestan 
province has been used. The result of the SPI drought 
index with an annual scale is shown in figure (1). 
According to the results of the SPI drought index, it can be 


—e— SPI 


So 

in 
= 
it 
= 


i=) 


SPI index 


S 
My 
—_—) 
j=) 
lee) 


year 


said that the Khorram River in the year 2012 was in 
2016 and 2018, 
respectively was wet and very wet. In other years, it has a 


moderate drought; in the years 


steady and near to normal trend. This indicates that the 
Khorram River has always been exposed to drought. 


1.76 


Figure 1. Result of drought in the years 2008 to 2018 in the study area 


3.2. WATER QUALITY INDEX RESULTS 


To calculate the WQI water quality index used from water 
quality data obtained from the Lorestan province regional 
water company. The results of the water quality index with 
an annual scale are shown in table (6). According to the 
results, it can be said that the water quality of the Khorram 
River has been in favorable condition according to the 


WQI index in the period of 2008 to 2018. This means that 
the water of the Khorram River was suitable for drinking, 
agriculture, and industrial uses. Table 7 shows the results 
of correlation coefficients of quantitative and qualitative 
parameters of Khorram River with influencing factors. 


Table 6. Values and description of WQI water quality index during the years 2008 to 2018 


Year WQI Index 
2008 82.23 
2009 75.38 
2010 71.13 
2011 74.75 
2012 74.46 
2013 71.76 
2014 80.13 
2015 75.27 
2016 73.59 
2017 70.60 
2018 72.12 


Descriptive 
equivalent 


good 
good 
good 
good 
good 
good 
good 
good 
good 
good 
good 


Table 7. Results of correlation coefficients of quantitative and qualitative parameters of Khorram River with influential 


Index SPI Precipitation TH PH 
WQi 
Quality 0.002 0.006 0.286 0.056 
Runoff 0.36 0.209 0.089 0.070 


factors 


Cl 


0.355 


0.14 
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Ca Mg SO, HCO; TDS Na 
0.267 0.295 0.076 0.170 0.672 0.151 
0.103 0.006 0.022 0.070 0.001 0.015 
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The result presented in table (7) shows that the highest and 
lowest correlation in quantitative parameters, respectively, 
is runoff-drought and runoff-total dissolved solids in water 
(TDS). Also, the highest and lowest correlation in quality 
parameters, respectively, is quality-total dissolved solids 


in water (TDS) and quality precipitation. The results of the 
correlation table do not show the communication of the 
quantities well, so used the quantities communication 
solution is based on coherence and cross wavelet. 


3.3. RESULTS OF IDENTIFYING THE RELATIONSHIP OF TIME SERIES BY COHERENCE AND 


CROSS WAVELET 


In this study, in order to determine the effect of climatic 
and human factors on changes in the quantity and quality 
of water in the Khorram river, we used four parameters 
precipitation, runoff, drought(SPI), and quality(WQD. 
Precipitation and drought are representative of climatic 
factors, and the runoff parameter due to increased human 
harvest from surface water resources to meet the needs of 
drinking, agriculture, and industry reduce the water flow is 
considered a representative human factor. Coherence and 
cross wavelet transformation have been used to evaluate 
the effect of the mentioned parameters on the quantity and 
quality of Khorram river water. Coherence and cross 
wavelet estimate the rate of interplay and phase delay of 
two-time series relative to each other. This transformation 
indicates in what period and with what phase delay the 
time series are related. For this purpose, the time series, 
quality (WQI) — precipitation, quality (WQI)-drought 
(SPI), runoff-precipitation, and runoff-drought (SPI) in 


pair enter the coherence and cross wavelet algorithm that 
is programmed in Matlab software, and are measured the 
rate of impact and interactions between these time series. 
In the results of coherence and cross wavelet, the areas 
marked with bold lines and dense arrows are the areas 
where the coherence and cross wavelet have shown local 
behavior. This concept means that there is a significant 
correlation between the corresponding fluctuations of the 
two time series in the periodicity (coherence band). The 
arrows to the right indicate that the two-time series are in 
the same phase, and the arrows to the left indicate that they 
are not in a phase. The up or down arrows indicate that a 
time series with a 90-degree angle leads to another. Figure 
3 shows the process of determining the impact of time 
series on the quantity and quality of Khorram river water 
by coherence and cross wavelet. The results of coherence 
and cross wavelet transformation for the time series are 


mentioned in figures 3, 4, and 5 are shown. 


Figure 2. The process of determining the impact of time series on the quantity and quality of Khorram river water by 


coherence and cross wavelet. 


A 
XWT: ‘Quality-Drought 


Period 


B 
XWT: 'Quality-Precipitation 


Period 


Figure 3. Results cross wavelet transformation for time series comparison, A: quality-drought and B: quality-precipitation. 
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WTC: ‘Quality-Precipitation 


Figure 4. Results coherence wavelet transformation for time series, A: quality-drought and B: quality-precipitation. 


Figure 4 shows the cross wavelet transformation of two 
quality-drought time series with a communication band 
range of approximately 16 to 32 months. On the other 
hand, according to figure 5, the wavelet coherence 
coefficient average obtained from this conversion is 
estimated to be 0.4. While the reciprocal wavelet 
transformation of two quality-precipitation time series 
with coherence wavelet band range from 16 to 32 months, 
the coherence wavelet coefficient is 0.2. The phase 
difference between the two quality-drought time series is 


A 
XWT: 'Runoff-Drought 


Sette 


20 40 60 80 100 120 


90 degrees according to the direction of the arrows, which 
are generally downward. Also, the phase difference of 
two-time series quality-precipitation according to the 
direction of the arrows, which are generally to the left, 
shows that the two-time series have opposed phases. Based 
on this interpretation, it can be said the relative effect of 
drought time series compared with precipitation time 
series on changes in water quality of Khorram River has 
been more from 2008 to 2018. 


B 
XWT: 'Runoff-Precipitation 


V4 


Figure 5. Results of cross wavelet transformation for time series comparison, A: Runoff-drought and B: Runoff- 
precipitation. 


A 


WTC: 'Runoff-Drought 


B 


WTC: ‘Runoff-Precipitation 


Figure 6. Results of coherence wavelet transformation for time series, A: runoff-drought and B: runoff-precipitation 
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Figure 5 shows the runoff-precipitation coherence wavelet 
transform with the communication band range of the two 
series for about 6 to 15 months. Another hand, the 
coherence wavelet average coefficient obtained from this 
transformation is estimated to be 0.6. In contrast, the 
coherence wavelet transformation of the two-time series 
runoff-drought with a band range of 16 to 32 months has a 
coherence wavelet coefficient of 0.1. Also, the phase 
difference between the two-time series runoff- 


precipitation and runoff-drought according to the direction 
of the arrow, which is generally to the right, shows that the 
two-time series are in the same phase. Based on this 
interpretation, it can be said the relative impact of 
precipitation time series compared with drought time 
series on changes in Khorram river runoff has been more 
from the years 2008 to 2018. Table (8) shows the results 
coherence wavelet average coefficient for hydrological 
time series. 


Table 8. Wavelet coherence coefficient average between Khorram river hydrological time series 


Time series Runoff-Drought 


Wavelet coherence coefficient average 0.1 


According to the results of table (8), it can be concluded 
that the low coherence wavelet coefficient between runoff- 
drought shows the effect of farmers’ uncontrolled 
withdrawals from Khorram river water and the effect of 
human factors in reducing the water runoff in Khorram 
River. On the other hand, the high coherence wavelet 
coefficient between runoff-precipitation shows the effect 
of climatic factors on the changes in water runoff in the 
Khorram River. Also, the high coherence wavelet 
coefficient of the quality-drought wave indicates that due 
to farmers' uncontrolled withdrawals from Khorram river 
water, the concentration of Khorram river water has also 
increased as a result. In a study, Yousefzadeh et al. (2013) 
studied the water quality of Khorram River using the water 
quality index (NSFWQI) and its zoning with GIS. Their 
results showed that the water quality of the Khorram River 


wavelet coherence coefficient and correlation 


Runoff-Precipitation Quality-Drought Quality-Precipitation 


0.6 0.4 0.2 


had the necessary standard for different uses [29]. It is 
worth mentioning that the most important parameter that 
is effective for accurately recognizing time series 
fluctuations by these two criteria (coherence and cross 
wavelet) is data continuity. Data continuity means that in 
the process of changing the time series of hydrological 
data, there is no discontinuity, so in order to maintain 
continuity during the time series of unrecorded data, there 
are no hydrological data during the time series 
fluctuations. On the other hand, coherence and cross 
wavelet criterion only deal with time series, so physical 
parameters such as soil permeability coefficients, type of 
vegetation in the area, etc., will not have any effect on the 
results of these two criteria. Figure 7 shows the coherence 
wavelet coefficient diagram and the correlation coefficient 
between time series. 


06 Green=Correlation Coefficient 
Orange=Wavelet Coherence Coefficient 


Figure 7. Diagram wavelet coherence coefficient average and correlation coefficient between time series 
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4. CONCLUSION 


Identifying the most important factors in lowering the 
quantity and quality of water resources is very important 
for future planning. In this research, first, drought analysis 
according to the SPI drought index based on precipitation 
data the 
organization has been done. The finding showed that in the 


from Lorestan province meteorological 
year 2012, the situation of the Khorram river was 
moderately drought in 2016 and 2018, the situation of the 
Khorram river respectively was moderately wet and very 
wet, and in other years it was near normal. Then the water 
runoff of the Khorram River and the effects of drought on 
it were studied. Study findings showed that the water 
runoff of the Khorram River has always been upward and 
downward and has not been stable in any year. It should be 
noted that this issue is one of the significant points and 
needs attention. The correlation between drought and 
runoff Khorram River with correlation coefficient is the 
scale of 0.36. in the next step, the water quality of the 
Khorram river was analyzed according to the global index 
(WQI), which the research findings showed the waters 
entering the Khorram river according to the standard of the 
world health organization (WHO) is in good, desirable 
conditions and suitable for various uses such as 
agriculture, drinking, industry and etc. using different 
methods and criteria to determine and classify the factors 
affecting the reduction of runoff and quality of aquifers is 
an effective step towards managing surface waters 
resources. In this regard, this research has tried to 


determine the relationship between different parameters 
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